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CN \ Abstract 

We deal with two following classes of equilibrium stochastic dynamics of infinite par- 
' ticle systems in continuum: hopping particles (also called Kawasaki dynamics), i.e., a 

dynamics where each particle randomly hops over the space, and birth-and-death pro- 
cess in continuum (or Glauber dynamics), i.e., a dynamics where there is no motion of 
particles, but rather particles die, or are born at random. We prove that a wide class 
of Glauber dynamics can be derived as a scaling limit of Kawasaki dynamics. More pre- 
cisely, we prove the convergence of respective generators on a set of cylinder functions, in 
Oh ' the L 2 -norm with respect to the invariant measure of the processes. The latter measure 

is supposed to be a Gibbs measure corresponding to a potential of pair interaction, in the 
low activity-high temperature regime. Our result generalizes that of [Finkclshtcin D.L. 
et al., to appear in Random Oper. Stochastic Equations], which was proved for a special 
Glauber (Kawasaki, respectively) dynamics. 



C3 



> ■ MSC: 60K35, 60J75, 60J80, 82C21, 82C22 

' Keywords: Birth-and-death process; Continuous system; Gibbs measure; Hopping particles; 

■ Scaling limit 

<N 

1 Introduction 

O 

This paper deals with two classes of equilibrium stochastic dynamics of infinite particle sys- 
tems in continuum. Let T denote the space of all locally finite subsets of R . Such a space 
is called the configuration space (of an infinite particle system in continuum) . Elements of T 
are called configurations and each point of a configuration represents position of a particle. 

One can naturally define a cr-algebra on T, and then a probability measure on T represents 
a random system of particles. A probability measure on T is often called a point process (see 
e.g. [9]). Configuration spaces and point processes are important tools of classical statistical 
mechanics of continuous systems. A central class of point processes which is studied there is 
the class of Gibbs measures. Typically one deals with Gibbs measures which correspond to 
a potential of pair interaction. 

An equilibrium stochastic dynamics in continuum is a Markov process on T which has a 
point process (typically a Gibbs measure) [i as its invariant measure. One can distinguish 
three main classes of stochastic dynamics: 

• diffusion processes, i.e., dynamics where each particle continuously moves in the space, 
see e.g. H El UM M M EU ; 
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• birth-and-death processes in continuum (Glauber dynamics), i.e., dynamics where there 
is no motion of particles, but rather particles disappear (die) or appear (are born) at 
random, see e.g. ^©01^113131131123]; 

• hopping particles (Kawasaki dynamics), i.e., dynamics where each particle randomly 
hops over the space [13j . 

For a deep understanding of these dynamics, it is important to see how they are related 
to each other. For example, in the recent paper [10], it was shown that a typical diffu- 
sion dynamics can be derived through a diffusive scaling limit of a corresponding Kawasaki 
dynamics. 

In [3], it was proved that a special Glauber dynamics can be derived through a scaling 
limit of Kawasaki dynamics. Furthermore, [3J conjectured that such a result holds, in fact, 
for a wide class of birth-and-death dynamics (dynamics of hopping particles, respectively), 
which are indexed by a parameter s € [0, 1]. (Note that the result of [3] corrrespods to the 
choice of parameter s = 0.) 

The aim of this work is to show that the conjecture of [3] is indeed true, at least for 
parameters s £ [0,1/2]. (In the case where s € (1/2,1], one needs to put additional, quite 
restrictive assumptions on the potential of pair interaction, and we will not treat this case 
in the present paper.) Thus, we show that the result of [3] is not a property of just one 
special Kawasaki (Glauber, respectively) dynamics, but rather represents a property which 
is common for many dynamics. 

More specifically, we fix a class of cylinder functions on T, and prove that on this class 
of functions, the corresponding generators converge in the L 2 (r,/i) -space. Here, /j, is a Gibbs 
measure in the low activity-high temperature regime, /U being invariant measure for all the 
processes under consideration. If one additionally knows that the class of cylinder func- 
tions is a core for the limiting generator, then our result implies weak convergence of finite- 
dimensional distributions of the corresponding processes. Unfortunately, apart from a very 
special case [11] , no result about a core for these generators is yet available. 

The paper is organized as follows. In Section 2, we briefly discuss Gibbs measures in the 
low activity-high temperature regime, and the corresponding correlation and Ursell functions. 
In Section 3, we describe classes of birth-and-death processes and of dynamics of hopping 
particles. In Section 4, we formulate and prove the result about convergence of the generators. 

The authors acknowledge numerous useful discussions with Dmirti Finkelshtein and Yuri 
Kondratiev. 

2 Gibbs measures in the low activity-high temperature regime 

The configuration space over M. d , d 6 N, is defined by 

r := {7 C M. d : [7 n A| < 00 for each compact A C M d }, 

where | • | denotes the cardinality of a set. One can identify any 7 G T with the positive Radon 
measure Ylx£-y £ x e -A4(M d ), where e x is the Dirac measure with mass at x, ^ £ % ■ — zero 
measure, and M(M. d ) stands for the set of all positive Radon measures on the Borel c-algebra 
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B(R d ). The space V can be endowed with the relative topology as a subset of the space A4(R d ) 
with the vague topology, i.e., the weakest topology on T with respect to which all maps 

r3 7 -></, 7 ):= / f(x) 7 (dx)=^2f(x), feC (R d ), 
J** x&1 

are continuous. Here, Co(R rf ) is the space of all continuous real-valued functions on R d with 
compact support. We will denote by B(T) the Borel cr-algebra on T. 

Let \i be a probability measure on (T,B(T)). Assume that, for each n £ N, there exists 
a non-negative, measurable symmetric function kj? on (M d ) Tl such that, for any measurable 
symmetric function /( n ) : (R rf ) n — > [0,+oo], 

/ ^ /" (n) (xi,...,x n )/i(d 7 ) = — / /W(a;i,...,a? fl )fcW(xi,...,x n )da;i---ds n . 

{xi,...,a; n }C7 v ' 

The functions A;^ are called correlation functions of the measure [i. If there exists a 
constant £ > such that 

V(xi, . . . , x n ) G (R d )« : fcW(xi, ...,x n )< r, (1) 

(n) 

then we say that the correlation functions k]i satisfy the Ruelle bound. 

The following lemma gives a characterization of the correlation functions in terms of the 
Laplace transform of a given point process, see e.g. [10] 



Lemma 1 Let fi be a probability measure on (T,B(T)) which satisfies the Ruelle bound ([I]). 
Let f : M. d — > R be a measurable function which is bounded outside a compact set A C M. d and 
such that e-F — 1 G L 1 (R d ,dx). Then for fi-a.a. 7 G F, (|/|, 7 ) < 00 and 



) . . . ( e /W - i)fc(») ( Xlj . . . , Xn ) dxi • • • cte n . 



Remark 1 Note that if / : R rf — > R is bounded outside a compact set A C R rf and if, 
furthermore, / is bounded from above on the whole R d , then the condition £ L 1 (M d , dx) 

is equivalent to / G L 1 (A c ,dx). 

Via a recursion formula, one can transform the correlation functions k^ into the Ursell 
functions u u and vice versa, see e.g. |22j . Their relation is given by 

kn{v) = ^u a (n 1 )---u tl (r] j ), r]£T o ,ri^0, (2) 

where 

r := {7 G r : | 7 | < 00}, 

for any rj = {x\, . . . , x n } E T 

ka{i) := A^ n) (xi,...,x n ), Up(rf) := u^(xi, . . . , x n ), 
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and the summation in ([2]) is over all partitions of the set r] into nonempty mutually disjoint 
subsets rji, . . . ,rjj C rj such that 771 U- • -Ur/j = rj, j G N. Note that if the correlation functions 
(k^)^ =1 are translation invariant, i.e., for each a G R rf 

k^{ Xl ,...,x n ) = kf\x l + a,...,x n + a), ( Xl ,...,x n ) G (R d ) n , 

then so are the Ursell functions (u^)^ =1 . 

A pair potential is a Borel-measurable function (ft : R rf — > R U {+00} such that 4>{—x) = 
(ft(x) G R for all x G R d \ {0}. For 7 G V and x G R d \ 7, we define a relative energy of 
interaction between a particle at x and the configuration 7 as follows: 



£(x, 7 ) := 



^0(x-y), if ^2\(/)(x-y)\ < +oo, 
+ 00, otherwise. 



A probability measure fi on (T,B(T)) is called a (grand canonical) Gibbs measure corre- 
sponding to the pair potential (ft an d activity z > if it satisfies the Georgii-Nguyen-Zessin 
identity ([T71 Theorem 2]): 

n(dj) / 7(dx)F(7, x) = / ^(dj) / zdxexpf— i?(x,7)]F(7 U x,x) (3) 

for any measurable function F : TxR^ — > [0; +00]. Here and below, for simplicity of notations, 
we just write x instead of {x}. We denote the set of all such measures [i by Q(z, (ft). 

As a straightforward corollary of the Georgii-Nguyen-Zessin identity ([3]), we get the 
following equality: 



H(df) I j(dxi) / y(dx 2 )F{i,xx,x 2 ) 

jR d JR d 

H(dj) / zdxi / zdx 2 ex.p[-E(xi,y) - E(x 2 ,y) - (ft(xi - x 2 )] 

JR d JR d 

x F(7 U {xi,x 2 },xi,x 2 ) 
+ / /J.(dj) / z dx exp [— E(x, 7)] F(j U x, x, x) (4) 



for any measurable function F : T x R rf x R rf — > [0, +00]. 
Let us formulate conditions on the pair potential (ft. 
(S) (Stability) There exists B > such that, for any 7 G To, 

^ 4>{x -y)> -B\~j\. 

{x,y}Cj 

In particular, condition (S) implies that (ft{x) > —2B, x G R rf . 
(LA-HT) (Low activity-high temperature regime) We have: 



/ \e~^ -l\zdx<(2e 1+2B )-\ 

jR d 
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where B is as in (S). 

The following classical theorem is due to Ruelle [21] 122] ■ 

Theorem 1 Assume that (S) and (LAHT) are satisfied. Then there exists fi G Q(z,(f>) which 
has the following properties: 

a) /x has correlation functions (A^)^Li, which are translation invariant and satisfy the 
Ruelle bound ([I]); 

b) For each n > 2, we have u^\o, •) G L 1 (M. d ^ n ^ 1 \ dx\ ■ ■ ■ dx n -\), where 
u^(0, •,-,...,•) is considered as a function of n — 1 variables. 

In what follows, we will assume that (S) and (LAHT) are satisfied, and we will keep the 
measure \x from Theorem [1] fixed. 

3 Equilibrium birth-and-death (Glauber) dynamics and hop- 
ping particles' 
(Kawasaki) dynamics 

In what follows, we will additionally assume that (f> is bounded outside some ball in M. d . Note 
that then (see e.g. [13]) 

E(x,j) = ^2<j)(x-y), 

2/67 

for dx //(c?7)-a.a. x G M. d and 7 G T and 

e(x, 7 \x)= Y, <P(x-y), 

y£f\x 

for fi-a.a. 7 G T and all x G 7. 

We fix a parameter s G [0, 1/2]. We introduce the set J-Cb(Co(M. d ), T) of all functions of 
the form 

T3 1 ^F( 1 )=g((f 1 , 1 ),...,(f N , 7 )), 

where N G N, /1, . . . , f N G C (M d ), and 5 G C^R^), where C^M^) denotes the set of all 
continuous bounded functions on M . . For each function F : T — > M, 7 G T, and x,y G 
we denote 

(D~F)(7) == F(7\x)-F( 7 ), 
p-+F)( 7 ) :=F(7\xUy)-F( 7 ). 

We fix a bounded function a : M rf — > [0, +00) such that a(— x) = a(x), cc G M d , and 
a G L 1 (M a! , dx). We define bilinear forms 

£ G (F,G) = [ M (d 7 ) / 7(^)exp[ Sj B(x,7\x)]( J D-F)( 7 )( J D-G)(7), 
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£k(F, G) = \! n{di) I j(dx) [ dy a{x - y) 



x exp[s£(x, 7 \ x) - (1 - s)E(y, 7 \ x)\{D~+ F)^)^ G)^) , 

where F,G<E J r C b (Co(M. d ),T). As we will see below, Sq corresponds to a Glauber dynamics 
and £k corresponds to a Kawasaki dynamics. 
The next theorem follows from |13j . 

Theorem 2 i) The bilinear forms (£ G ,J r C b (C (R d ),r)) and {S K , FC b (C (R d ) ,T)) are clos- 
able on L 2 (T, /j.) and their closures are denoted by (£g, D{£g)) an d {£ki D(£k)), respectively. 

ii) Denote by (Hg,D(Hg)) and (Hk,D(Hk)) the generators of (£g,D(£g)) and 
(£ K ,D(£ K )), respectively. Then J r C b (C (R d ),T) C D{H G ) n D{H K ), and for any F € 

FC b (c (R d ),r), 

(H G F)( 7 ) = - [ 7 (dx)exp[ Sj E(x, 7 \ a ;)]( J D-F)( 7 ) 

-[ zdxexp[(a-l)E(x,i)](D+F)(rt), (5) 

JR d 

{H K F){ 1 ) = - f 7 (<fx) / dya(x-y)eMsE(x,^\x) + (s-l)E(y, 1 \x)}(D-+F)( 1 ). 

JR d JR d 

(6) 

hi) Let (j := G,K. There exists a conservative Hunt process 

M» = F», (F{) t > , (0f) t >o, (X«(t) t > , (P^) 7 er) 

on T (see e.g. [151 P- 92]) which is properly associated with (£j, D(fjj)) ; i.e., for all {^-versions 
of) F 6 L 2 (r,/i) and all t > the function 



r9 7 H+(p»F)( 7 ):= / F(X«(t))dP« 



is an fjj- quasi- continuous version of ex.~p[—tH$]F . M" is np to /i- equivalence unique (cf. 
Chap. IV, Sect. 6]). In particular, M" /ias /x as invariant measure. 

Remark 2 In Theorem can be taken canonical, i.e., fi" is the set -D([0, +oo), T) of all 

cddldg functions u : [0, +oo) — » T (i.e., cj is right continuous on [0, +oo) and has left limits 
on (0, +oo)); X"(i)(u;) = u;(i), i > 0, u G fi"; (Ff) t > , together with F", is the corresponding 
minimum completed admissible family (cf. [5, Section 4.1]); ©|, i > 0, are the corresponding 
natural time shifts. 

It follows from (0) that Hq is (at least heuristically) the generator of a birth-and-death 
process, in which the factor exp[si?(x, 7 \ x)\ describes the rate at which particle x of the 
configuration 7 dies, whereas the factor exp[(s — l)E(x, 7 )] describes the rate at which, given 
a configuration 7 , a new particle is born at x. We see that particles tend to die in high energy 
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regions, i.e., if E(x,'y \ x) is high, and they tend to be born al low energy regions, i.e., if 
E(x, 7) is low. 

Next, by ©, Hk is (again at least heuristically) the generator of a hopping particle 
dynamics, in which the factor 



describes the rate at which a particle x of configuration 7 hops to y. We see that this rate is 
high if the relative energy of interaction between x and the rest of the configuration, 7 \ x, 
is high, whereas the relative energy of interaction between y and 7 \ x is low, i.e., particles 
tend to hop from high energy regions to low energy regions. 

4 Scaling limit 

In this section, we will show that the birth-and-death dynamics considered in Section 3 can 
be treated as a limiting dynamics of hopping particles. In other words, we will perform 
a scaling of Kawasaki dynamics which will lead to the Glauber dynamics. We will only 
discuss this convergence at the level of convergence of the generators on an appropriate set of 
cylinder functions. In fact, such a convergence implies weak convergence of finite-dimensional 
distributions of corresponding equilibrium processes if additionally the set of test functions 
forms a core for the limiting generator. However, in the general case, no core of this generator 
is yet known and this is an open, important problem, which we hope to return to in our 
future research. Our results will hold for all s E [0, 1/2] (see Section 3). They will generalize 
Theorem 4.1 in [3], which was proved in the special case s = 0, and confirm the conjecture 
formulated in Section 6 of that paper. 

So, let us consider the following scaling of the Kawasaki dynamics (for a fixed s G [0, 1/2]). 
Recall that, for each bounded function a : R d — > R such that a(x) > 0, a G L 1 (M. d , dx), and 
a(— x) = a(x) for all x G R rf , we have constructed the corresponding generator of the Kawasaki 
dynamics. We now fix an arbitrary e > and define a function a £ : M. d — > R by 



By the properties of the function a, we evidently have that the function a £ is also bounded, 
satisfies a E (x) > 0, for all x G R d , a £ G L 1 (R d ,cix), and a £ {— x) = a £ {x) for all x G R d . 
Hence, we can construct the Kawasaki generator which corresponds to the function a £ . It is 
convineant for us to denote this generator by (H £ , D(H £ )). We will also denote the generator 
of the Glauber dynamics by (Hq, D(Hq)). We first need the following lemma, whose proof is 
completely analogous to the proof of Lemma 4.1 in [3]. 

Lemma 2 For any e > and any if G C (R d ), the function F(j) := e^' 7 > belongs to D(H £ ) 
and the action of H £ on F is given by the right hand side of formula © for e > (with a 
replaced by a £ ), respectively by the right hand side of ([5]) for e = 0. 



exp[s£(x, 7 \ x) + (s - l)E(y, 7 \ x)] 



a £ (x) = e d a(ex) 



Note that 




a(x) dx. 
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Remark 3 For each e > 0, denote by (£ £ , D(£ £ )) the Dirichlet form with the generator 
(H £ , D(H £ )). It can be easily proved that the set { exp[(<y9, •)] : tp G Co(K rf )} is dense in the 
Hilbert space D(£ £ ) equipped with inner product (F, G)^^) '■= £(F, G) + (F, G) L 2rp ifJ/ y 

We have 

(D~ y + F)( 1 ) = F( 1 \xUy)-F( 1 ) 

= -F( 7 ) + F( 7 \ x) - F( 7 \ x) + F( 7 \ x U y) 
= (D~F)( 1 ) + (D+F)( 1 \x). 

So, we may rewrite the action of H £ for e > as follows: 

H E :=H+ + H-, 

where 

{HJF){ 1 ) = - f 7 (da;)(Z»-F)(7) / dna £ (x-y)exp[sL;(£, 7 \x) + (s-l)L%, 7 \x)] 

jR d JR d 

and 

(H+F)( 7 ) = -/ 7 (dx) / dy a e (a?- J /)exp[*^(a;,7\x) + (8-l)£(y,7\x)](D+F)(7\x). 
We can also rewrite 

#0 := #o + -^o ' 

where 

(F -F)( 7 ) = - / 7 (dx)exp[ S £?(x,7\x)](Z)-F)(7) 

JR d 

(F +F)( 7 ) = - / Z dxexp[-(1 - S )£(x, 7 )](I?+F)(7). 

JR d 



and 



Theorem 3 Lei s G [0,1/2] be fixed. Assume that the pair potential <p and activity z > 
satisfy conditions (S) and (LA-HT). Assume that <f> is bounded outside some compact set in 
M. d . Assume also that 

4>{x) — ► as \x\ — » oo. (7) 

Lei /x 6e £/ie Gi66s measure from Q(z,<p) as in Theorem^ Assume that the function a is 
chosen so that 

-l 



/ a(x)dx = ( 




/R d V 





T/ten, /or any <p G C (lR d ), 

fl^e^ -► H+e^ in L 2 (r,/i) as e -> 0, 

so i/iai 

ff e e^>'> -► fToe^'") in L 2 (r, /i) ase^ 0. 



M (d 7 ) • (8) 
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Remark 4 In fact, condition ([7]) can be omitted, and instead one can use the fact that (ft 
is an integrable function outside a compact set in IR rf (compare with [3]). However, in any 
reasonable application, the potential (ft does satisfy condition (|7|). 

Remark 5 Note that the integral on the right hand side of ([8]) is well defined and finite due 
to Lemma [U see also Remark [T] 

Proof. We first need the following lemma, which generalizes Lemma 4.2 in [3]. 

Lemma 3 Let a function ift : M. d — > M. be such that e^ — 1 is bounded and integrable. Suppose 
that A > 0, B > 0, xi, x 2 , j/i, y 2 € M. d and x\^y\. Then 



exp 



AE\^ + x 2 ,i 



exp 



fj,(dj) /^exp 



m£7 



K d l) / exp[(^,7)]M(^7) 



as e —* 0. 

Proof. By Lemma [TJ 



exp 



/i(d7) 



- AE((xi/e) +x 2 ,7) -BE({ yi /e) +y 2 ,j) + (^,7> 
°° 1 /" 

1 + E "I / (exp[-^(- - - - y(e)) + V(0] - l)®(t*i, . - - , n n ) 
£zj n! J (R d)n 

x k^\ui, . . . ,u n )dui ■ ■ -du n , (9) 



where x(e) := (xi/e) + x 2 , y{e) := (yi/e) + y 2 . 

Using the Ruelle bound, semi-boundedness of (ft from below and the integrability of (ft 
outside a compact set, we conclude from the dominated convergence theorem that, in order 
to find the limit of the right hand side of ([9]) as e — ► 0, it suffices to find the limit of each 
term 



: 



(exp[— A(ft(- — x(e)) 



x kjfi (u\ , . . . , u n ) du\ ■ ■ ■ du n 

E ( n 

ni+n2+n 3 =n x 

x kp>(ut, . . . ,u n ) dui ■ ■ ■ du n , 



y(e)) + i>(-)}-l)® n ( Ul ,...,u n ) 



{f?^®f^T®ft*)(ui,...,u n ) 



(10) 



where 



h,e( u ) '■= ( ex Ph^(«)] - l)exp[-Acft(u - x(e)) - B<ft{u - y(s))], 
f2,e(u) ■= (exp[-A/>(u - x(e))] - 1) exp[-B(ft(u - y(e))], 
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f 3;£ (u) := exp[-S0(n - y(e))] - 1, u£t. 
Using definition of Ursell functions, we see that 

(AT ® / 2 T ® / 3 T)K, • • • » «n)^ n) • • • , «n) • • ■ 

£ / (AT ® AT ® /s? 3 )^ ■ ■ ■ , MMvi) ■ ■ ■ Mm) ■ ■ ■ du n , 



where the summation is over all partitions {771, . . . , rjj} of 77 = {u\, . . . ,u n }. We now have to 
distinguish the three following cases. 

Case 1: Each element rfc of the partition is either a subset of {u±, . . . , u ni }, or a subset of 
}, or a subset of {u ni+n2+ i, . . . , u n }. Set 

u i = Ui - x(e), i = rti + 1, . . . , n 2 , 

u'i = Ui — y(e), i = n 2 + 1, . . . , n. 

Then using the translation invariance of the Ursell functions, we get that the corresponding 
term is equal to 

(ffe 1 ®9ie 2 ® gfe 3 )(ui, ... ,11^11^) ■■■ U^rjj) dll! ■■■ du n , (11) 

J (R d ) n 

where 

92,e(u) ■= {exp[-A(j)(u)] - 1) exp{-Bcf)(u + ((xi - y\)/e) + x 2 - J/2)], 
g 3>E (u) := exp[-B<t>(u)] - 1, u G M. 

Note that x\ — y\ ^ and so for any fixed u (and x 2 , y 2 ), we have 

I it + ((xi - yi)/e) + x 2 - y 2 \ — > +00 as e — > 0. 

By ([7]) and the dominated convergence theorem, we therefore have that (jll j) converges to 

(exp[-</'(-)] - If" ® (exp[-A0(-)] - 1) 



/ 

J(R d ) 



(exp[-B </)(•)] - l)®" 3 (iti, . . . , MM 7 !!) ■ ■ ■ Mm) dui--- du n . 



Case 2: There is an element of the partition which has non-empty intersections with both 
sets {ui, . . . ,u ni } and {u ni+1 , . . . ,u n }. 

Using Theorem [IJ we have that, for each n € N, 

jjjn+i) € L 1 ((R d ) n ,dx 1 ---dx n ), 

where 

U^ +1 \ Xl ,...,x n ) :=4" +1 )(xi,...,x n ,0), (xi,...,x„) G (M d ) n 
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Consider the integral 

/ (exp[—0(ui)] - l)u^(ui, ...,u k ) dui ■■■duk 

J(R d ) k 

(exp[-V>(ui)] - l)uff\0,u 2 - U!,U3 -ui,...,u k - ui) dui ■ ■ ■ du k , 



(R d ) k 

where we used translation invariance of Ursell function. By changing variables u[ = u±, u' 2 = 
U2 — u\, . . . , u' k = u k — ui, we continue as follows: 

= / (exp[-V>K)] - (0, u' 2 , u' 3 , . . . , u' k ) dui • • • H 

= / (e"^ (Ml) - x / U( k \ u2 ,u 3 ,...,u k )du 2 ---du k . 

jR d J(R d ) k - 1 



Note also that 

\x(e)\ — > +00 and |y(e)| — > +00 as e — > 0, 
and hence, for each fixed u €.M. d 

exp[-Acp(u - x{e))} - 1 -» 0, exp[-B^(u - y(e))] - 1 

as e — > 0. From here, using the dominated convergence theorem, we conclude that 

J(R d )™ 

x Upim) ■ ■ ■ UpiVj) dui--- du n -> 

as e — > 0. 

Case 3: Case 2 is not satisfied, but there is an element 77/ of the partition which has 
non-empty intersections with both sets {u ni+ i, . . . , u ni+n2 }, and {u ni+n2+ i, . . . , u n }. 

Shift all the variables entering r]i by x(e). Now, since exp [-A<f>] - 1 G dx), analo- 

gously to case 2, the term converges to zero as e — > 0. 

Thus, again using the definition of the Ursell functions, we get, for each n G N, 

y f n ^ 

x / (exp[<0(-)] -l) 0ni (^i ) --->"n 1 )A:J t ni) (ui,...,M ni )dui---du ni 

J(R d ) n l 

x / (exp[-A/>(-)] - l)® n2 {u ni+1 , . . . ,u ni+n2 )k^ 2 \u ni+ i, . . . ,u ni+n2 ) 

J(R d ) n 2 



x du ni +\ • ■ ■ dii rai + n2 



x / (exp[-B0(-)] - l) 0n3 (u ni+n2+ i, . . . ,u n )k^ 3 \u ni+n2+1 , . . . ,u r , 

J(R d ) n 3 

x du ni+ri2+ i ■ ■ ■ du n . 
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Therefore, the right hand side of (|9j) converges to 

f oo . „ 

l + V- / (eM-M-)}-l) m (ui,...,u n )k^(u 1 ,...,u n )du 1 ---du 

1 + / (exp[-S^(-)] -l)® n (ui,...,u n )kfr\ui,...,u n )d Ul ---du n ) 

(l + Y,-, f (eKp[iP(-)]-l) m (u 1 ,...,u n )kP(u 1 ,...,u n )du 1 ---du n ) 



/r 6XP 




ji(d 7 ) y" exp 




M (d 7 ) y 








L «G7 





as e — > 0, which proves the lemma. □ 

Now we are in position to prove the theorem. We fix any <p E Co(R ) and denote ^(7) := 
e \<psv _ it suffices to prove that 

j{H£Ff{j)ii{drt) - J (H±F) 2 (j)pL(dj) as e - 0, (12) 



o 2 



(H?F)(>y)(H*F)( 1 )»(d'y) - / (H^Fy^idj) as 6 -> 0. (13) 
r Jr 

Now, 

j^Ff^nm 

( - [ j(dx)exp[sE(x^\x)}(D-F)( 7 )) M (d 7 ) 
ridl)( - [ 7(dx)exp[sE(x, 1 \x)]{e^ ) -^ x) - e { ^) 

\ JR d , 

M(d7) f - / 7(*=) exp[s^(x, 7 \ x)]e^' 7) (e" 1 ^ - 1)^) 

^(d7)e^' 7> / 7 (dxi) / 7(dx 2 )exp[s£(xi,7\xi)]exp[s£ , (x2,7\x 2 )] 
x (e-^ 1 ) - l)(e~^ 2) - 1) 

fi(drf) f zdxe- E(x ^e {2 ^ fUx) exp[2s£(x, 7 )]( e -^ ) - l) 2 

JR d 

+ /" M(d 7 ) f zd Xl [ zdx 2 exp[-E(x 1 ,7)-E(x2,7)-0(xi-x 2 )]e {2 ^ UxiUX2) 
Jr JR d JR d 

x exp[s J B(xi,7Ux 2 )]exp[sF(x 2 ,7Uxi)](e-^ (:Cl) - l){e-^ X2) - 1) 
/ 2 dx(l - e^f I fjt(dry) exp[(2s - l)E(x, 7) + {2ip, 7)] 
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r 



+ / zdxi zdx 2 e^ xl \l-e^)e^ X2 \l-e^) 

JR d JR d 

x exp[(2s - l)0(xi - x 2 )] j M(d7)exp[(s - l)E{x ul ) + (s - l)E(x 2 ,-y) + (2<A7>] 

(14) 

and 

r (flfi0 2 (7)M(*y) 

zdx x [ zdx 2 (e^ xi) - l)(e^ X2) - 1) / ^7) exp[(s - l)E( Xl ,^) 
jR d Jr 

+ {s-l)E{x 2 , 1 ) + {2^ 1 )\. (15) 
Completely analogously to (fT4"j) we have 

(H-F) 2 ( 7 )fi(d 7 ) 

zdx dy% I dy 2 (e ip ^ - l) 2 a e (x - y\)a e (x - y 2 ) \ /J,(dj) 
jR d jR d Jr 

exp[(2s - l)E(x,-y) + {s- l)E(y in ) + {s - l)E(y 2 ,j) + (2<p, 7)] 

+ / zdx! f zdx 2 [ d Vl f dy 2 e^ Xl) (e^ Xl) -l)e^ X2 \e^ X2) -1) 

JR d jR d JR d JR d 

x a £ (xi - yi)a e (x 2 - y 2 ) exp[(2s - l)<j){xi - x 2 ) + (s - l)<j>(yi - x 2 ) 
+ (s - l)<f>(xi - y 2 )\ J n(d-y) exp[(s - l)E( Xl , 7) + (s - l)E{x 2 ,i) 
+ {s- l)E( yi ,j) + (s- l)E(y 2 ,i) + (V, 7)]- (16) 
Let us make the change of variables 

y'i=£(yi-x), y' 2 = e{y 2 -x) 

in the first integral, and 

y[ = e(y 1 -x 1 ), y' 2 = e{y 2 -x 2 ) 

in the second integral. Then omitting the primes in the notations of variables, we continue 
(fT5|) as follows: 



zdx / d yi / dy 2 (e^ - l) 2 a{ yi )a{y 2 ) / ^7) exp[(2s - l)E(x, 7 ) 
jR d JR d Jr 

+ (s- l)E(( yi /e) + x, 7) + (s - l)E((y 2 /e) + x,j) + (2<p, 7}] 

+ / zdxi ! zdx 2 [ dyi I dy 2 e^ie^ -l)e v{x2) (e^ X2) -1) 

JR d JR d JR d JR d 

x a(yx)a{y 2 ) exp[(2s - l)<j>(x\ - x 2 ) + (s - l)(f){{yi/e) +x x - x 2 ) 
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+ (s - l)0((jft/e) + x 2 - xi) y n(dj) exp[(s - l)£(xi, 7) + (s - l)E(x 2 , 7) 

+ (a - l)£(( yi /e) + xi, 7) + (s - lMfo/e) + x 2 , 7) + (2^, 7}]. (17) 

Next, 

V+F) 2 (7)/x(cZ 7 ) 

= / H(<h)( - / 7(<**0 / dya £ {x-y) 
Jr \ JR d JR d 

\ 2 

x exp[sE(x, 7 \ x) - (1 - s)E(y, 7 \ x)](F( 7 \ x U y) - F( 7 \ x)) 

= / m(^7)( / 7(dx) / dya £ {x-y) 
Jr \ JR d JR d 

x exp[s^(x, 7 \ x) - (1 - s)E(y,j\ x )]( e ^'^ + ^ - e ^M*>) 
= I »{d 7 )( [ 7 (dxi) / 7 (cix 2 )e^\ a;i >e^\ :E2 > / / dy 2 

Jr \JR d JR d JR d JR d 

x a £ (x 1 -yi)a £ (x 2 - y 2 )exp[sE(x 1 ,^\x 1 ) - (1 - s)E(y ± , 7 \ xi)] 
x exp[s£(z 2 , 7 \ x 2 ) - (1 - s)E(y 2 , 7 \ x 2 )](e^ - l)(e^ - 1) 

= / M(<*7) / 7(^)e <2v ' 7 ^ > / dj/i / dy 2 

x a e (x - yi)a £ (x - j/ 2 )exp[s£(x,7\ x) - (1 - s)E(y 1 , 7 \ x)] 
x exp[sS(x, 7 \x) - (1 - s)E(y 2 ,~f \x)](e^ - l)(e^ - 1) 

+ / nidi) I 7(dxi) / ( 7 \x 1 )(dx 2 )e^ :El >e^^ a;2 > / dyi / dy 2 

Jr JR d JR d JR d JR d 

x a e (xi -yi)a £ (x 2 - y 2 )exp[sE(x 1 ,^\x 1 ) - (1 - s)E(y ± , 7 \ xi)] 
x exp[ Sj E(x 2 ,7\x 2 ) - (1 -s)E(y 2 ,7 \x 2 )](e^ - l)(e^ - 1) 

= / M(<*7) / ^dxexp[- J B(x,7)]e< 2 ^ / dyi / dy 2 
Jr JR d JR d JR d 

x e M a(e(x - yi))a(e(x - y 2 )) exp[sE(x,-/) - (1 - s).E(?/i, 7)] 
x exp[s£(x, 7 ) - (1 - s)E(y 2l 7)] (e v(i/l) - l)(e<^ - 1) 

+ / /i((i7) / zdxi / zcfe 2 / dyi / dy 2 

Jr JR d JR d JR d JR d 

x exp[- J E(xi, 7 ) - £(x 2 , 7 ) - </>( Xl - x 2 )]e^ UX2) e { ^ Uxi) 
x a e (xi -yi)a E (x 2 - y 2 ) expfs^^i, 7 U x 2 ) - (1 - s)E(y 1 , 7 U x 2 )] 
x exp[s J E(x 2 ,7Ux 1 ) - (1 - s)%,7Uii)](e v(!/l) - l)(e^ 2) - 1) 
= : I + II. (18) 
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In the first integral in (|18p let us make the change of variables 

y'i = £(vi-x), y'2 = e {v2-x). 

Then omitting the primes in the notations of variables, we continue I as follows: 

1= f dy x I dy 2 a(yx)a( y y 2 )(e^ yi ^ +x ^ - l)( e ¥>((w/e)+s) _ i) 

JR d JR d 

x / fj,(drf) / zdx exp[(2s - l)E(x,i) + (s - l)E{{ yx /e) + x,j) 
+ {s-l)E((y 2 /e)+x n ) + (2<p, 1 )]. 

Let us take 

x' = x + (yi/e), 
then omitting the primes in the notations of variables, we get: 

1= / d Vl f dy 2 a( yi )a(y 2 ) [ z dx(e^ - l){ e ^+^-yi)/e)) _ i) a ( yi ) a ( y2 ) 

JR d JR d JR d 

x J fi(drf) exp[(2s - l)E(x - ( yi /e), 7) 
+ (s - l)E{x, 7) + (s - l)E(x + ((y 2 - yi)/e), 7) + (2<p, 7)]. 
In the second integral in (|18p , let us make the change of variables 

x' 1 =e(x 1 -y 1 ), x' 2 = e(x 2 -y 2 ). 
Then omitting the primes, we have: 

11=/ zdx x ! zdx 2 [ d Vl I dy 2 e' fiiixi/£)+yi) e' fi{{x2/£)+y2) (e viyi) - l)(e^ fe) - 

JW. d JR d JR d JR d 

x a(xx)a(x 2 ) exp[(2s - 1)</>(((xi - x 2 )/e) + yx - 2/2) + (s - l)0((a:i/e) + yi - 
+ (s - l)0(yi - y 2 - (x 2 /e))] jT M (d 7 ) exp[(s - l)£( yi , 7) + (a - 7) 
+ (a - l)^((x!/e) + y 1>7 ) + (s - l)E((x 2 /e) + 3/2,7) + (2ip, 7)]. 
Using (HU) and ([20]), we get 

| r (F+F) 2 ( 7 )/i(d 7 ) 

= / zdx/ dyi / dy 2 (e^ - i)( e ^+((s«-s/i)A) _ l)o(yi)o( W ) 

JM d JR d JR d 

x jT /i(d 7 ) exp[(2s - l)£(x - (yi/s), 7) 

+ (s - l)E(x, 7) + (s - l)E(x + ((y 2 - yi )/s), 7) + (2y>, 7}] 
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+ / zdxxf zdx 2 [ dm f dy 2 e^ xl ^ +yi ^e^ X2 ^ +y2 \e^ - l)(e^ - 1) 

JR d JR d JR d JR d 

x a(x 1 )a(x 2 ) exp[(2s - l)0(((xi - x 2 )/e) +y x - y 2 ) + (s - l)0((xi/e) + yi - y 2 ) 
+ (s - l)0(yi - 2/2 - (afe/e)) jf m(^7) exp[(s - l)£(yi,7) + (s - l)£(y 2 , 7) 
+ (a - l)^((ari/e) + j/i, 7) + (a - + 2/2, 7) + <2y, 7>]- (21) 

Completely analogously, we get 

^(if -F)(7)(ff-F)( 7 ) M (d 7 ) 

= / z / (e<^) - l) 2 a(2/) / M (d 7 ) exp[(2s - l)E(x, 7) + (s - l)E{{y/e) + x, 7) + <2y>, 7)] 
jR d Ju d Jr 

+ [ zdx x \ zdx 2 [ dy{e^ Xl) - l)(e^ 2) - l)e^ Xl) a{y) 

JR d JR d JR d 

x exp[(2s - l)4>(xi - x 2 ) + (s - l)0((y/e) + »i - x 2 )] 

x y M (d 7 ) exp[(s - l)£(xi,7) + (« - 1)^(^2,7) + (« - l)E((y/e) + a*, 7) + (2</>,7>], 

(22) 

and 

= [ zdxi [ zdx 2 [ dya(x 2 )(e^ x ^ - 1) 

JR d JR d JR d 

x ( e v(i/) _ i) jf M (d 7 ) e (2^,7> exp[(s _ l)£(( X2 /e) + y, 7) + (s - l)£(y, 7) 

+ (a - l)E(xi, 7) + (s - l)H(x 2 /e) +y- Xl )]. (23) 

Using the Ruelle bound and Lemma [U we conclude that the integral over T in the right 
hand side of equalities (fTTj) . (f2"Tj) - (f23"j) . are bounded by a constant, which is indepeandent of 
e. Therefore, by the dominated convergence theorem, to find the limit of (|17p . (|21|) - (|23p as 
e — > it suffices to find the point-wise limit of the functions appearing before the integral 
over r, as well as the limit of the integrals over T for fixed x (x\ and x 2 respectively), y% and 

To find the latter limits, we use Lemms [3l Then, using ([8]), we see that (|17p and (|22p 
converge to ()14p . whereas ()21[) and (|23p converge to (|15p . Therefore, (|12p and Q13p hold. □ 
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